Susceptibility and Percolation in 2D Random Field Ising Magnets 
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The ground state structure of the two-dimensional random field Ising magnet is studied using 
exact numerical calculations. First we show that the ferromagnetism, which exists for small system 
sizes, vanishes with a large excitation at a random field strength dependent length scale. This 
break-up length scale Lb scales exponentially with the squared random field, exp(j4/A^). By adding 
an external field H we then study the susceptibility in the ground state. If L > L;,, domains melt 
continuously and the magnetization has a smooth behavior, independent of system size, and the 
susceptibility decays as L~^. We define a random field strength dependent critical external field 
value ±i/c(A), for the up and down spins to form a percolation type of spanning cluster. The 
percolation transition is in the standard short-range correlated percolation universality class. The 
mass of the spanning cluster increases with decreasing A and the critical external field approaches 
zero for vanishing random field strength, implying the critical field scaling (for Gaussian disorder) 
He ^ {A - Ac)^ where Ac = 1.65 ± 0.05 and 5 = 2.05 ± 0.10. Below Ac the systems should 
percolate even when H — 0. This implies that even for H — above Lb the domains can be fractal 
at low random fields, such that the largest domain spans the system at low random field strength 
values and its mass has the fractal dimension of standard percolation Df = 91/48. The structure 
of the spanning clusters is studied by defining red clusters, in analogy to the "red sites" of ordinary 
site-percolation. The size of red clusters defines an extra length scale, independent of L. 
05.50.+q, 75.60. Ch, 75.50.Lk, 64.6O.Ak 
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I. INTRODUCTION 



The question of the importance of quenched random 
field (RF) disorder in ferromagnets traces back to the 
primary paper by Imry and Ma |^J^ from mid-seventies. 
They argued using energy minimization for an excitation 
to the ground state, that the randomness in the fields as- 
signed to spins changes the lower critical dimension from 
the pure case with d/ = 1 to d/ = 2. After that a number 
of field-theoretical calculations suggested that the ran- 
domness increases d/ with two to be d/ = 3. Finally 
came rigorous proofs first by Bricmont and Kupiainen |^] 
in '87 that there is a ferromagnetic phase in the three- 
dimensional (3D) random field Ising model (RFIM) and 
in '89 by Aizenman and Wehr Q that there is no fer- 
romagnetic phase in 2D RFIM. Thus it was established 
that the lower critical dimension is two. This means that 
the ground state is a paramagnet, but the problem as to 
how to describe the structure of the (ground state of) 2D 
RFIM still persists. Some recent work concerns the scal- 
ing of the correlation lengths ^ and there is a suggestion 
of a ferromagnetic phase, but with a magnetization that 
is in the thermodynamic limit below unity . The point 
is that due to the (relevant) disorder there are no easy ar- 
guments that would indicate, say, how the paramagnetic 
ground state should be characterized. This is different 
from the thermal Ising case, which is quite trivial in ID. 

In two dimensional Ising magnets, in the presence 
of quenched random fields, the problem of determining 
the ground state (GS) becomes more difficult. Finding 
the true ground state with any standard Monte Carlo 
method is problematic due to the complex energy land- 



scape. Even with the exact ground state methods, as 
the one used in this paper, the thermodynamic limit is 
difficult to reach, since the finite size effects are strong. 
In the typical case of square lattices the only way not 
to have a massive domain, which would scale with the 
Euclidean dimension of the system, is to have enough in- 
terpenetrating domains of both spin orientations. How- 
ever, with decreasing strength of the randomness the 
ferromagnetic coupling constants between spins start to 
matter, the domains become "thicker" , and thus one en- 
ters an apparent ferromagnetic regime, and the paramag- 
netic (PM) phase is encountered only at very large length 
scales. Should there be large clusters with a fractal (non- 
Euclidian) mass scaling they nevertheless can contribute 
to the physics in spite of the fact that the total fraction 
of spins can be negligible in the thermodynamics limit. 
Thus such clusters may even be measurable in experi- 
ments or be related to the dynamical behavior in non- 
equilibrium conditions. Therefore it is of interested to 
study the structure of the large(st) clusters in the ground 
state, since it is not simply paramagnetic like in normal 
Ising magnets above Tc. The true ground state structure 
gives also some insight into the physics at T > 0, since 
the overlap between the GS and the corresponding finite- 
T state is close to unity for T small, in contrast to the 
thermal chaos in spin glasses [Q. 

In this paper we want to shed some light on the char- 
acter of the ground states of 2D RFIM. We have done 
extensive exact ground state calculations in order to char- 
acterize how the ferromagnetic (FM) order vanishes with 
increasing system size. We have also studied the effect of 
the application of an external field, that is the suscepti- 
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bility of the 2D RFIM. Allowing for a non-zero external 
field makes it possible to investigate a percolation -type of 
critical phenomenon for the largest clusters. We propose 
a phase diagram in the disorder strength and external 
field plane for the percolation behavior. The presence of 
clusters of the size of the system, i.e. percolation type 
of order brings another correlation length in the systems 
and thus makes the decay of ferromagnetic order more 
complicated than at first sight. 

The Hamiltonian of the random field Ising model is 



n 



(1) 



where J > (in this paper we use J = 1 for numerical 
calculations) is the coupling constant between nearest- 
neighbor spins Si and Sj. We use here square-lattice. H 
is a constant external field, which if non-zero is assigned 
to all of the spins, and hi is the random field, acting on 
each spin Si. We consider mainly a Gaussian distribution 
for the random field values 
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with the disorder strength given by A, the standard de- 
viation of the distribution, though in some cases the bi- 
modal distribution, 



A) + -Sih, 



A) 



(3) 



is used, too. The results presented below should not be 
too dependent on the actual P{h), in any case. 

To find the ground state structure of the RFIM means 
that the Hamiltonian (|l|) is minimized, in which case the 
positive ferromagnetic coupling constants prefer to have 
all the spins aligned to the same direction. On the other 
hand the random field contribution is to have the spins 
to be parallel with the local field, and thus has a para- 
magnetic effect. This competition of ferromagnetic and 
paramagnetic effects leads to a complicated energy land- 
scape and the finding the GS becomes a global optimiza- 
tion problem. An interesting side of the RFIM is that 
it has an experimental realization as a diluted antifer- 
romagnet in a field (DAFF). By gauge-transforming the 
Hamiltonian of DAFF 



TC — — J SiSjEiej 



B ^ eiSi. 



(4) 



where the coupling constants J < 0, is the occupation 
probability of a spin Si, and B is now a constant ex- 
ternal field, one gets the Hamiltonian of RFIM (|^) with 
H ^0 §-[131. The ferromagnetic order in the RFIM cor- 
responds to antiferromagnetic order in the DAFF, natu- 
rally. 

As background, it is of interest to review a few basic 
results. Imry and Ma used a domain-wall argument to 



show that the lower critical dimension d/ = 2 [Q . In or- 
der to have a domain there is an energy cost of 0{L'^^^) 
from the domain wall. On the other hand the system 
gains energy by fiipping the domain from the fluctuations 
of random fields, which interpreted as a typical fluctua- 
tion means that the gain is 0{L'^/^). Thus, whenever 
d/2 > d — 1, i.e. d < 2, it is energetically favorable for 
the system to break into domains. However, in this pa- 
per we will point out, as it has been shown in ID [ pT| , 
that the 0{L'^/'^) scahng can be used only in relation to 
the sum of the random fields in the "first excitation" , but 
not to the droplet field energy when the GS consists of 
domains of different length scales. 

Grinstein and Ma derived from the continuum in- 
terface Hamiltonian that the roughness of the domain 
wall (DW) in RFIM scales as w ^ ^2/32^(5-d)/3^ ^i-^^^i-^ 
is consistent with di = 2. Later Fisher used the 
functional renormalization group (FRG) to obtain the 
roughness exponent C = (5 — d)/3 and argued that due 
to the existence of many metastable states the pertur- 
bative RG calculations and dimensional reduction fail. 



Another, microscopic calculation by Binder |14| opti- 
mized the domain- wall energy in two-dimensions. The 
net result is a total energy gain from random fields, 
AU = — (A^/J)LlnL due to domain wall decorations, 
which implies that the domain wall energy U ~ 2JL+AU 
vanishes on a minimal length scale 



Lb^exp[A{J/A)% 



(5) 



where A is a constant of order unity. For L > the 
expectation is that the system spontaneously breaks up 
into domains. Similarly the energy of a domain with a 
constant external field H becomes 

U/J = 8L + 2{H/J)L^ - {8/A){A/J)^LliiL. (6) 

Setting U/J = and assuming that the critical length 
scale Lb^H scales as Lb, without the field, i.e. Lb.n ^ Lb, 
the critical external field becomes 



i7,/J= (4M)(A/J)2exphA(J/A)2 



(7) 



Note, that in this case the first two terms in (||) assume 
that the domain is compact. 

These results imply, together with the notion that 
the ground state is paramagnetic, that the magneti- 
zation should not as such display any "universal" fea- 
tures. The results of this paper show that the magneti- 
zation is not dependent on the system size and has a 
smooth behavior of ni ~ f[H/ exp(— 6.5/A)] and the 
susceptibility vanishes with the system size as x ~ 
L-2exp(7.3/A)g[i//exp(-6.5/A)], where g{y ~ 0) ~ 
const, and g{y ±00) ^ exp(— 0.2|?/|). These imply 
that there is a length scale, related to the rate at which 
clusters "melt" when H is changed from zero. 
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The presence of such a length scale is qualitatively sim- 
ilar to the one discovered in the context of the percolation 
transition. It turns out so that when the external field is 
varied, the universality class is that of the ordinary short- 
range correlated percolation universality class. The ex- 
ternal field threshold for spanning He with respect to de- 
creasing random field strength approaches zero external 
field limit from the site-percolation limit of infinite ran- 
dom field strength value, suggesting a behavior for Gaus- 
sian disorder of 7?^ ~ (A - Ac f , where Ac = 1.65 ± 0.05 
and S = 2.05±0.10. Below this value the lattice effects of 
site percolation are washed out, and there is yet another 
length scale that characterizes the percolation clusters, 
the size of the "red clusters" defined below in analogy to 
the usual red or cutting sites in percolation. Now a whole 
cluster is reversed due to the forced reversal of a "seed" 
cluster, when the sample is optimized again. The length 
scale is however finite, indicating that the global opti- 
mization of the ground state creates only finite spin-spin 
correlations as is the case in ID as well. 

This paper is organized so that it starts by introduc- 
ing in Section |l^ the exact ground state calculation tech- 
nique. In Section [II the breaking up of ferromagnetic 
order is discussed, based on a nucleation of droplets - 
picture which follows from a level crossing between a FM 
ground state and one with a large droplet. The relevant 
Lfc scaling (||) is derived from extreme statistics. Above 
the break-up length scale the domains have a complex 
structure that is briefly discussed. The effect of an exter- 
nal field, in the case the system size is above L^, is studied 
in Section IV for Gaussian disorder. The percolation as- 
pects of the 2D RFIM are studied in detail in Section 0. 
The phase diagram for the percolation behavior as func- 
tions of the external field and the random field strength 
is sketched and the properties of the transition are dis- 
cussed. The zero-external field percolation probability 
is studied in Section VI. In the section also the struc- 



ture of spanning clusters is studied using the so called 
red clusters whose scaling and properties are discussed. 
The paper is finished with conclusions in Section VII. 



II. NUMERICAL METHOD 

For the numerical calculations the Hamiltonian is 
transformed to a random flow graph with two extra sites: 
the source and the sink. The positive field values hi 
correspond flow capacities cu connected to the sink (t) 
from a spin Si, similarly the negative fields with Cis are 
connected to the source (s), and the coupling constants 
Jij = 2cij between the spins correspond fiow capacities 
Cij = Cji from a site Si to its neighboring one Sj . In 
the case the external field is applied, only the local sum 
of fields, H + hi, is added to a spin towards the direction 
it is positive. The graph-theoretical combinatorial op- 
timization algorithms, namely maximum-flow minimum- 
cut algorithms, enable us to find the bottleneck, which 



restricts the amount of the flow which is possible to get 
from the source to the sink through the capacities, of 
such a random graph. This bottleneck, path P which di- 
vides the system in two parts: sites connected to the sink 
and sites connected to the source, is the global minimum 
cut of the graph and the sum of the capacities belong- 
ing to the cut Cij equals the maximum fiow, and is 
smaller than of any other path cutting the system. The 
value of the maximum flow gives the total minimum en- 
ergy of the system. The maximum flow algorithms are 
proven to give the exact minimum cut of all the random 
graphs, in which the capacities are positive and with a 
single source and sink p6[ |. In physical situations this 
means the systems are without local frustration. The 
algorithm was actually used for the first time in this con- 
text by Ogielski , who showed that the 3D RFIM has 
a ferromagnetic phase. The best known maximum flow 
method is by Ford and Fulkerson and called the aug- 
menting path method . We have used a more sophis- 
ticated method called push-and-relabel by Goldberg and 
Tarjan | |T9[ |, which we have optimized for our purposes. It 
scales almost linearly, 0(n}'^), with the number of spins 
and gives the ground state in about minute for a million 
of spins in a workstation. 

When we have added an external field in the systems 
our system sizes are restricted to = 175^, for H small 
but nonzero due to the range of integer variables (for 
numerical reasons we use a discrete representation of 
real fields). When the high precision for field values is 
not needed the computations extend up to system sizes 
— 1000^. We have used periodic boundary conditions 
in all of the cases. Also the percolation is tested in the 
periodical or cylindrical way, i.e., a cluster has to meet 
itself when crossing a boundary in order to span a system. 



When the red clusters are studied. Section VIC, we 



have applied a technique, which allows us to take ad- 
vantage of the so called residual graph After the 
original ground state is searched a perturbation is ap- 
plied. This means that e.g., a spin is forced to be to the 
other direction with a large opposite field value. Then 
the ground state is searched again. This time all the flow 
need not to be constructed from scratch, but instead one 
can utilize the final situation of the first ground state 
search (the residual graph). Only the extra amount of 
fiow, needed since the capacity of the large opposite field 
value is added, has to be forced through the system to 
the sink. One has to subtract also the flow from the orig- 
inal field value (retrace it back to the source). It is thus 
convenient to reverse only the fields which originally were 
negative. For the positive field values one would have to 
study a mirror copy of the system (hi — > —hi). Thus we 
have analyzed the red clusters only from the spanning 
clusters of down spins, which does not disturb the statis- 
tics, since the spin directions are symmetrical. The use 
of the residual graph reduces the time to calculate the 
next ground state considerably, although approximately 
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a half of the spanning clusters have to be neglected. No- 
tice that since the ground state energy is a linear function 
of the capacity of the saturated bonds (or the field of the 
spins aligned along the local field) one can compute the 
"break-point" field hb, at which a change takes place from 
the original ground state to the new one. We have not 
paid attention to this, however, due to our main interest 
in the geometry of the red clusters. One interesting ad- 
ditional question would be what is the smallest hb and 
its disorder-averaged distribution. 



III. DESTRUCTION OF FERROMAGNETIC 
ORDER 

In this section we will derive the scaling for the break- 
up length scale Lt, Eq. (^), from extreme statistics (as 
done in the paper by Emig and Nattermann [|T]). and 
confirm it with exact ground state calculations. We also 
discuss the ensuing domain structure qualitatively. 

If one picks a (compact) subregion of area a of a fer- 
romagnetic 2D RF system the energy is drawn from a 
Gaussian distribution 

where the variance tr^ — aA^ is due to the fluctuations 
of random fields, and (E) ~ a. For a system of size 
we have Na ^ ways of making such a subregion. The 
probability that a subregion has the lowest energy E is 
given by 

LNAE)^NaViE){l-CiiE)}''--\ (9) 

where CiiE) = f^^P{e)de |2|. The distribution 
Lfq^{E) is in fact a Gumbel distribution The av- 

erage value of the lowest energies is given by 

/oo 
ELMAE)dE, (10) 
-oo 

which can not be solved analytically. The typical value 
of the lowest energy follows from an extreme scaling esti- 
mate. The factor inside the curly brackets in (||) is close 
to unity if Ci becomes small enough (for similar applica- 
tions, see |§|-§^). Thus 

aNaV{{Eo)) ~ 1 (11) 

which yields, 

{Eo) ^ {E) - <j{\n{Na)V/^ (12) 
i.e. the energy gain from the fluctuations is 

{E,)^a{HNa)V/'. (13) 



A FM system would tend to take advantage of such 
large favorable energy fluctuations by reversing a domain, 
which requires breaking bonds. This is assumed to have 
a cost of 

Eb - Ja(''-i)/f (14) 

Equating Eqs. ( p^ and ( pH ) yields the estimate of the 
parameter values at which the flrst "Imry-Ma" domain 
occurs, 

v^A{ln(iVa)}^''^ Ja^'^-i)/'^. (15) 

It can be easily understood that the most preferable do- 
main is the one, which maximizes the area and minimizes 
the bonds to be broken, which gives a ~ I? 12. Fig. |l| il- 
lustrates this, as we increase (with a fixed random field 
configuration and system size) the strength of the ran- 
domness or decrease the ferromagnetic couplings until 
the first domain appears. It turns out to be so that the 
droplet is of the order of the system size. This kind of nu- 
cleation with a critical size is reminiscent of a first order 
transition, and is related to a level-crossing, when either 
the random field strength or the system size is varied, 
similarly to random elastic manifolds, when an extra pe- 
riodic potential |^ or a constant external field ||2^ is 
applied. By substituting a ^ and Na ^ to (|15|), 
we get for the length scale 

L ~ exp [AiJ/Af] , (16) 

which is in fact Lb as in (^). This result, (|l6|), is surpris- 
ing in the sense, that the extreme statistics calculation 
for the formation of a domain leads to the exactly same 
scaling as the optimization of a domain wall energy on 
successive scales in Binder's argument. 

Due to the extensive size of the first domain-like ex- 
citation the destruction of the ferromagnetism resembles 
a first order transition. The magnetization for a cer- 
tain disorder strength and system size would be averaged 
over systems, in which the excitation has and has not 
been formed yet, with |m| ~ and |?7i| ~ 1, respectively. 
Hence we define a simpler measure for the break up of FM 
order: the probability of finding a purely ferromagnetic 
system, Ppm{L, A), i.e., for a fixed random field strength 
and system size we calculate the probability over several 
realizations that magnetization |r7i| = 1 |^^. If the tran- 
sition to the PM state would be continuous, this would 
not make much sense, since already small fluctuations 
would cause Ppm{L,A) ~ 0. However, due to the first 
order behavior, and to the fact that the smallest energy 
needed to flip a domain causes the excitation to be large, 
PpM is a good measure and has a smooth behavior. We 
have checked that \m\ vs. Pfm does not depend on L. 

We have derived the break-up length scale Lb by vary- 
ing the random field strength A from the probability of 
finding a pure ferromagnetic system as PpAfiLb, A) = 
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0.5. The data is shown in Fig. || for Gaussian and bi- 
modal disorder (in both cases J = 1), and the expo- 
nential scahng for Li, vs. inverse random field strength 
squared is clearly seen. The prefactors are A = 2.1 ± 0.2 
and 1.9 ± 0.2 for Gaussian and bimodal disorder, respec- 
tively. To check that the probability Pfm{L, A) < 1 is 
not due to so called stiff spins, i.e., single spins for which 
hi > AJ , we next derive an extreme statistics formula 
for their existence. The probability of finding hi > AJ 
P{hi > 4 J) = erfc(4/A). The extreme statistics argu- 
ment, NP{h, > 4J) « 1, with N ^ L'^, gives 

Li « [erfc(4/A)]-i/^ (17) 

For Gaussian A ~ 0.8862 for which Lj, = 100 in Fig. | 
Li ~ 400 from Eq. (|l7|). Li also grows much faster 
than Lb for decreasing A which is both easy to see from 
Eq. (^^ and to check numerically. For A ~ 0.6670 for 
which Lb ~ 800 Li becomes as huge as 22300. To con- 
firm further that the origin for the break-up is a large 
domain, one can extend the argument to small domains. 
The length scale L2 at which one is able to find a cluster 
of two neighboring spins flipped, i.e., N/2 L'^P{hi + h2 > 
6 J) 1, where hi < 4 J, 6J — hi < h2 < hi, becomes 
even greater than Li. These small clusters are present 
in large system sizes, but do not play a role in the prob- 
ability of first excitation, since the energy minimization 
prefers extensive domains. It is amusing to note that 
the critical droplet size reminds of critical nucleation in 
ordinary first order phase transitions. It is also worth 
pointing out that the reasoning for stiff spins does not 
work for the bimodal distribution (since it is bounded), 
and indeed we observe as expected a similar Lb scaling 
for both the Gaussian and bimodal disorders. 

When a system size is well above the break-up length 
scale the Imry-Ma argument is no longer applicable to the 
structure. In Fig. |3| we depict two systems with a large 
Gaussian random field strength value A ~ 3.0 and with 
a smaller one A = 1.2 for a system size L^ — 100^. One 
can see that a system breaks to smaller and smaller do- 
mains inside each other from the case seen in Fig. |l|. The 
feature of having clusters in different scales is familiar 
from the percolation problem 1 29 1 . In fact in both of the 
examples in Fig. |^ there is a domain, which spans the sys- 
tem in vertical direction, drawn in gray. For the stronger 
random field value one can see also smaller domains of 
different sizes. However, the width of the spanning clus- 
ter is greater than in a standard site or bond occupation 
percolation problem. Later in subsection VI B we discuss 



the scaling properties of the largest clusters in the ground 
state, above Lb- 



IV. MAGNETIZATION AND SUSCEPTIBILITY 
WITH AN EXTERNAL FIELD 

In Fig. ^ we show what happens in system, with a 
system size well above Lb when an external field H is 



applied. Now the clusters "melt" smoothly when the 
external field strength is increased and such first order 
type of phenomenon is not seen as when a first Imry-Ma 
droplet appears in the zero-field case. The magnetization 
has a continuous behavior, see Fig. |^(a), where we have 
the magnetization with respect to the external field for 
several Gaussian disorder strength values. All the mag- 
netization values for different system sizes lie exactly on 
top of each other, when L > Lb, and as long as the statis- 
tics is good. 

For smaller system sizes L < Lb one could study 
"avalanche" -like behavior, i.e. the sizes of the areas that 
get turned with the magnetic field (see [Q). However, 
these are due to the first order break-up, defined by (|^) 
and one should bear in mind that such behavior does 
not exist in the thermodynamic limit, L — > 00, when the 
system sizes are above the break-up length scale. For 
L > Lb our results indicate that the size distribution of 
the ffipped regions as H is swept is not that interesting. 

In order to find the scaling between the external 
field and the random field strength we have taken from 
Fig. |5|(a) the crossing points of magnetization curves with 
a fixed magnetization values at external fields Hm for dif- 
ferent random field strength values A. The external field 
Hm scales exponentially with respect to the random field 
strength. 



~ exp(-6.5/A), 



(18) 



see Fig. ^(b). This is also an evidence of non-existence 
of a critical point in A, in which case there should be a 
power-law behavior if the transition was continuous, thus 
no PM to FM transition is seen. The data-collapse us- 
ing the scaling (|l^) is shown in Fig. ^(c) confirming the 
prediction of the scaling. The magnetization has a linear 
behavior with respect to the external field for small field 
values H and exponential tails. The exponential behav- 
ior of Eq. ([1^) implies that there is a unique "melting 
rate" at which the cluster boundaries get eroded as H 
increases and that the process is otherwise similar for all 
A. We have no analytical argument for the melting rate 
[or the slope of the m{H)-cmve], and note that it is not 
seemingly, at least, related to Lb. 



We have also studied the susceptibility, x 



(m 



{m)'), with respect to the external field. In Fig. ||(a) the 
susceptibility is shown for a fixed random field strength 
A = 2.2 and varying system size. The data has been 
collapsed with the area of the systems, x/^~^- I^i 
Fig. H(b) we have data-collapsed the susceptibility ver- 
sus random field strength by scaling the external field 
with (^8|) as for magnetization and the susceptibility with 
X ~ exp(7.3/A). Again the exponential behavior is a sign 
of non-existence of any critical point, due to the lack of 
power-law divergence at any Ac. Although the shape of 
the data-collapse of the susceptibility looks almost Gaus- 
sian, it is actually not. It has a constant value for small 
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external field values H and exponential tails for large val- 
ues, as seen in Fig. ||(c). This results straightforwardly 
from the magnetization, since x = dm/dH. To summa- 
rize the behavior of the susceptibility, it is 



X^L"2exp(7.3/A).g(i//iI™), 
where is from Eq. (fsh and 



9iy) 



const, 
exp(-0.2|y| 



0, 

> ±00. 



(19) 



(20) 



Therefore the fluctuations of the magnetization are asso- 
ciated with yet another scale, which is almost but not 
quite an inverse of that related to the magnetization. 
From the suscpetibility one gets the magnetization cor- 
relation length which has an exponential dependence 
on the random field strength. It should be noted, finally, 
that we have here studied only the case with Gaussian 
disorder. With any other distribution we would expect 
that the prefactors in Eqs. (|l^), (p^, and ( pO| ) would 
change. 



V. PERCOLATION WITH AN EXTERNAL FIELD 

Motivated by Fig. H, where the domains resemble the 
percolation problem we next study the percolation be- 
havior in the 2D random field Ising magnets with Gaus- 
sian disorder. The usual bimodal distribution could be 
studied as well, but since it is susceptible to some anoma- 
lous features we concentrate on the Gaussian case which 
does not have these problems. The bimodal case suffers 
from the fact that the ground states arc highly degenerate 
at fractional field strength values. Thus there are some 
ambiguities in how define percolation clusters [^,|2| . 

When the random field strength is well above the cou- 
pling constant value, A 3> J, the percolation can be 
easily understood by considering it as an ordinary site- 
occupation problem. This means that, only the ran- 
dom field directions are important and the coupling con- 
stants may be neglected. The site-percolation occu- 
pation threshold probability for square-lattices is Pc — 
0.593 [|2^, i.e., well above one half. Applied to the strong 
random field strength case it means that there must be 
a finite external field in order to get a domain spanning 
the system. However, when the random field strength 
is decreased, the coupling constants start to contribute 
and in some cases systems span even without an external 
field, as in Fig. ||. Hence, we propose a phase-diagram 
Fig. 1^. There we can take the limit 1/ A = so that the 
ordinary site-percolation problem is encountered. This is 
true for distributions for which one can control the frac- 
tion of "stiff" spins (i.e. hi > 4 J) systematically. In the 
case of Gaussian disorder there will be of course, even 
for A very large a small fraction of "soft" spins where 
this criterion is not fulfilled. Thus the exact point that 



the percolation line approaches in the 1/ A ^ 0-limit will 
depend on the distribution, but we expect that the "pc" 
is different from one-half, that is He ^ 0. Notice that 
again the binary distribution presents a problem. 

When 1/A 00 the percolation threshold lines start 
to approach each other and the H = line. Now there 
arises two questions. The first one is, what kind of a 
transition is the percolation here? Is it like the ordinary 
short-range correlated percolation, suggested by the site- 
percolation analogy for the strong random field strength 
case, or are there extra correlations due to the global op- 
timization relevant here? Examples about similar cases 
can be found from Ref. |^^. The second question is, do 
the lines meet each other at finite Ac, i.e. does there 
exist a spanning cluster also when H — and A > 0? 
Our aim is to answer these two questions in this section, 
where we study the percolation problem in the vertical 
direction in the phase-diagram Fig. |^, and in the next 
section, where the horizontal direction, H — line, is 
considered. 

In Fig. ||(a) we have drawn the spanning probabili- 
ties of up spins Hup with respect to the external field 
H for several system sizes L, which are greater than Li,, 
for a fixed random field strength A = 2.6. The curves 
look rather similar to the standard percolation. When we 
take the crossing points Hc{L) of the spanning probabil- 
ity curves with fixed spanning probability values for each 
systems size L, we get an estimate for the critical external 
field He using finite size scaling, see Fig. ||(b). There we 
have attempted successfully to find the value for He using 
the standard short-range correlated 2D percolation cor- 
relation length exponent v = 4/3. Using the estimated 
He ~ 0.00186 for A = 2.6 we show a data-collapse of Hup 
versus [H — He)/L~^/^ in Fig. p(c), which confirms the 
estimates of He and v — 4/3 |29|. We get similar data- 
collapses for various other random field strength values 
A as well. In order to test further the universality class 
of the percolation transition studied here, we have also 
calculated the order parameter of the percolation, the 
probability of belonging to the up-spin spanning cluster 
Poo ■ Using the scaling analysis for the correlation length 



^perc ^ I LI He \ , 

and for the order parameter, when L < £,perc, 

Poo{H) ^ {H ~ Hef , 

we get the limiting behaviors, 

{H — HeY L < £,perc, 

T -P/^ T ^ f 

^ ^ ^ Sperc: 



Poo{H,L) 



(21) 
(22) 

(23) 



and thus the scaling behavior for the order parameter 
becomes 



Poo{H,L)^ L-P/-F 



{H-He)-" 
L 



(24) 
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We have done successful data-collapses, i.e. plotted 
the scaling function /, for various A using the stan- 
dard 2D short-range correlated percolation exponents 
(3 = 5/36 and v — 4/3, of which the case A — 3.0 with 
He = 0.040 is shown in Fig. ||. Note, that the values 
{H—Hc)/L^^'''^ > are not shown, since cut-offs appear, 
due to the fact that Poo is bounded in between [0, 1] and 
the non-scaled Poo — 1 values after scaling saturate at 
different levels depending on the system size. 

Hence, we conclude that the percolation transition for 
a fixed A versus the external field H is in the standard 2D 
short-range correlated percolation universality class . 
This is confirmed by the fractal dimension of the span- 
ning cluster, too, as discussed below. Also other expo- 
nents could be measured, as 7 for the average size (s) 
of the clusters, and a and t for the cluster size distri- 
bution. Note, however, that then the control parameter 
should be the external field H instead of the disorder 
strength A. See Ref. jsj] for an example of the clus- 
ter size distribution for a non-critical case (H — 0, but 
\Hc\ ^ 0). The other exponents should be measurable, 
too, like the fractal dimension of the backbone of the 
spanning cluster, the fractal dimension of the chemical 
distance, the hull exponent etc. In addition to the cor- 
rect control parameter also the break-up length scale, has 
to be considered, too. Notice that there is an slight con- 
tradiction hidden in the notion that the hull exponent 
could be measured. Namely, both the work of Ref. fs^ j 
and studies of domain walls enforced with appropriate 
boundary conditions give no evidence thereof. It seems 
likely that to recover the right exponent (4/3) one has to 
resort to studying the spanning cluster geometry itself at 
the critical point with L > L^. It is amusing to note that 
the standard 2D percolation hull exponent can be recov- 
ered in non-equilibrium simulations of 2D RFIM domain 
walls 35 1 . 

We have now shown, that there exists a line of critical 
external field values -ffc(A) for percolation. The corre- 
sponding correlation length ^perc diverges as (|2^) with 
a correlation length exponent v = 4/3. On the other 
hand it was shown in the previous section, that there is 
no critical external field value for the magnetic behavior, 
i.e., no PM to FM transition, and the magnetic corre- 
lation length bas an exponential dependence on A. 
The percolation correlation length ^perc may cause some 
confusion, when studying the PM structure of the GS, 
since it brings in another length scale. 

To answer the question how the percolation critical ex- 
ternal field i?c(A) behaves with respect to the random 
field, we have attempted a critical type of scaling using 
the calculated He for various A =2.0, 2.2, 2.6, 3.0, and 
4.5. For smaller A Lt becomes large and He approaches 
the vicinity of zero, being thus numerically difficult to 
define. We have been able to use the Ansatz behavior of 



the calculated A values versus the scaled critical ex- 
ternal field [iJc(A)]^/^ °^ and it gives the estimate for 
Ac = 1.65 ± 0.05. This indicates that the percolation 
probability lines for up and down spins meet at Ac — 1.65 
and for A below the critical Ac there is always in the sys- 
tems spanning of either of the spin directions, even for 
H = 0. Actually one should note, that the only way that 
neither of the spin directions span is to have a so called 
checker-board situation, which prevent both of the spin 
directions to have neighbors with the same spin orienta- 
tion. However, the another scenario with an exponential 
behavior for Hc{A) fits also reasonably well. This would 
suggest, that there is no finite Ac. Fig. |l^(b) shows a 
behavior of He ^ A'^ cxp{-13/ A'^ - 4). This can be 
compared with Eq. (|^) , where the break-up external field 
was derived. Notice that the derivation was for compact 
domains and the spanning clusters here are by default 
fractals. Besides that, the factor 13 in front of 1/A^ is 
much larger than A = 2.1 in the scaling form for Lf,. The 
difference implies that the Lc, at which length scale the 
spanning probability vanishes, scales as Le ^ L^. The 
is already exponentially large length scale for small A, so 
Le should be large enough that one can be below it in 
experiments, and thus a system can "apparently perco- 
late" 1^0. 



VI. PERCOLATION AT H = 

To understand how the percolation transition is seen 
when there is no external field and the random field 
stren gth is changed we study the phase diagram in Sec- 
tion VIA in the direction of the horizontal arrow in 
Fig. [7|. The str uctur e of th e span ning clusters are stud- 
ied in Sections VI B and in VI C with the help of the so 
called red clusters. 



He^iA- AeY 



(25) 



where S = 2.05 ± 0.10. In Fig. |q(a) we have plotted 



A. Spanning probability 

In Fig. |l^(a) we have plotted the probability for span- 
ning of either up or down spins H as a function of the 
Gaussian random field strength A. The probabilities are 
calculated up to A = 30, but only the interesting part 
of the plot is shown. There is a drop from H = 1 at A, 
which corresponds L = Lj, for each system size, to a value 
about H ~ 0.85. We have also calculated the H^p in this 
case and it is approximately one half of H. For the larger 
random field strength values the probabilities H decrease 
and the lines get steeper, when the system size increases. 
In order to see, if the spanning probabilities are converg- 
ing towards a step function at some threshold value, we 
have calculated the probabilities up to the system size 
= 1000^ and each point with 5000 realizations. 

For each system size L we have searched the cross- 
ing points Ac{L) of the spanning probability curves in 
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B The percolation cluster 



Fig. |ri](a) with fixed probabihty values H =0.1, 0.15, 0.2, 
0.25, 0.3, 0.4, and 0.5. Using finite size scahng for Ac{L) 
of the form Ac{L) ^ A^l + CiL-^^''){l + G^^^/"^) we 
have estimated Ac for each 11 value, see Fig.[ll|(b). There 
we have plotted the Ac(£) versus the scaled system size 
CiL-^/". One sees that the different threshold values for 
spanning probabilities 11 approach different critical ran- 
dom field strength values Ac. The threshold II's have 
been plotted with respect to Ac in Fig. |ll](c). In the 
ordinary percolation, this should be a step function, and 
the correlation length exponent v independent on the cri- 
terion n. However, here also l/v is dependent on the cri- 
terion n and varies with respect to Ac (11). We believe, 
that this surprising phenomenon is due to that we are ap- 
proaching the part in the phase diagram Fig. 0, where the 
percolation lines of up and down spins are getting close 
to each other. In terms of the two control parameters A 
and H one can think about the "percolation manifold" : 
it has a line of unstable fixed points i/c(A). Usually H 
is a good control parameter close to H^- Having A as 
a control parameter seems to have the problem that one 
moves almost parallel to the actual line Hc{A). 

When considering the percolation probability of up or 
down spins, it actually consists of probabilities of up- 
spin spanning n„p and down-spin spanning lidown as 
n = Hup + Udown - n^pHdo^.n, since they are corre- 
lated with each other. Assuming that H„p (and Hdown 
respectively) has a value about one-half at the critical 
line of percolation at the thermodynamic limit, we get 
n = 0.75. In standard percolation such a value is not 
actually universal (and we have not confirmed it) but de- 
pends on the boundary conditions, etc. [ ^Tf . However, 
whatever the values for n„p and Hdown are at the ther- 
modynamic limit, as long as they are below unity, H is 
below unity, too. This may be the reason, why there is 
an immediate drop in Fig. |ll](a) from H = 1 at A{Lb) 
for each system size, to a value about H ~ 0.85. If we 
approximate with a linear behavior the H versus Ac(H) 
in Fig. Pt|(c), the critical value estimated in the previous 
section Ac = 1.65 ± 0.05, when the percolation threshold 
He = for up-spin spanning, has a value about H ~ 0.7. 
Another interesting point in Fig. ^(c) is that the l/v 
is about 3/4 , when H versus AjTl) approaches zero. 
Thus the standard correlation length exponent would be 
reached far enough away from the area, where the perco- 
lation threshold lines for up and down spins touch each 
other. 

In order to test the break-up length scale type of scal- 
ing for the percolation behavior (|^) , we have taken from 
Fig. |l|(a) the estimated Ac(L) for various H values and 
plotted the system sizes in double-logarithm scale ver- 
sus the logarithm of the inverse of the critical Ac(£), see 
Fig. ^(d). The exponent, which is a = 2 in Lb scaling, 
Lb ~ exp(A/A"), is now dependent on H again. At least 
this does not solve the problem here, and the break-up 
length scale type of scaling can be ruled out. 



B. The percolation cluster 

In order to see if the thickness of the spanning cluster 
affects the scaling of the standard percolation we have 
measured the fractal dimension of the spanning cluster 
when H = 0. By now unsurprisingly, the standard two- 
dimensional short-range correlated percolation fractal di- 
mension Df = 91/48 fits very well in the data, as can be 
seen in Fig. |2|. The least-squares fit gives a value of 
Df — 1.90 ±0.01. We have measured also the sum of the 
random fields in the spanning cluster and found that the 
sum scales with the exponent Df — 91/48, too. This is 
in contrast to the Imry-Ma domain argument, where the 
sum is taken scale as L''/^. The prefactor for the scal- 
ing of the sum of the random fields approaches slowly 
zero with decreasing random field strength, opposite to 
the mass of the spanning cluster, which increases with 
decreasing A. 

Hence, the Imry-Ma argument defines only the first 
excitation, and is irrelevant when it comes to domains 
when the system has broken up to many clusters on dif- 
ferent length scales. Then the structure is due to a more 
complicated optimization. The domains are no longer 
compact and as noted above for large enough domains 
the domain-wall length should be characterized by the 
percolation hull exponent. 



C. Red clusters 

So far all the evidence points out to the direction that 
the percolation transition is exactly of the normal uni- 
versality class. To further investigate the nature of the 
clusters in the presence of the correlations from the GS 
optimization, we next look at the so-called red clusters. 
The structure of a standard percolation cluster can be 
characterized with the help of the "colored sites" picture 
in which one assesses the role of an element to the connec- 
tivity of the spanning cluster. This picture has been also 
called the links-nodes-blobs model with dead-ends [ p9| . 
The red sites, or links and nodes, are such that removing 
any single one breaks up the spanning cluster. 

To compare with the original ground state we investi- 
gate what happens if one inverts, by fixing the local field 
hi to a large value opposite to the spin orientation, any 
spin belonging to the spanning cluster. Then the new 
GS is found with this change to the original problem. 
The effect is illustrated in Fig. |^, the crucial difference to 
site-percolation is that now a whole sub-cluster can be re- 
versed. The spin drawn in yellow is the inverted, "seed" 
spin in the spanning cluster, and the spins painted red 
form the rest of the red cluster, which are flipped from the 
original ground state when the energy is minimized the 
second time. We do the investigation whether the origi- 
nal cluster retains its spanning property for each spin or 
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trial cluster in analogy with ordinary percolation. Those 
spins that lead to a destructive (cluster) flip, define then 
red clusters (RC) as all the spins that reversed simulta- 
neously. 

The finite size scaling of the number of the red clus- 
ters, (Nrc), is shown in the Fig. |l^(a) for different field 
values. (Njic) is in practice calculated as the number of 
the seed spins, which cause the breaking of the spanning 
cluster, since two different seed spins may both belong 
to the same red clusters without the red clusters being 
identical. The technique to find the red clusters was in- 
troduced in Section || and although it is efficient, only 
up to the system size 0(200^) can be studied, since each 
of the spins in spanning clusters has to be checked sep- 
arately, because one cannot know beforehand, whether 
it is critical or not - this is what we want to find out. 
For smaller field values the spanning cluster is "thicker" 
and the red clusters get larger. One can see from the 
Fig. |l|(a) that (Nrc) scales with L^/", where ~ 4/3 
as in ordinary percolation, for field values A < Ac, when 
L > Lf,. The amplitude is larger the smaller the field, 
as is the average mass of red clusters {Mjic)- {M^c) is 
independent of the s^tem size L and depends only on 
the field A, see Fig. ||(b). 

The other elements of the spanning cluster, dead-ends 
and blobs, could be generalized, too. Here blobs, which 
are multiple connected to the rest of the spanning cluster, 
are such that in order to break a spanning cluster, sev- 
eral seed spins are needed to flip simultaneously instead 
of a single one. Links, nodes, and blobs form together 
the backbone of the spanning cluster, and the rest of the 
mass of the cluster is in the dead-ends. The red cluster 
size scale defines the average smallest size of any element 
of the spanning cluster. 



VII. CONCLUSIONS 

In this paper we have studied the character of the 
ground state of the two-dimensional random field Ising 
magnet. We have shown that the break-up of the ferro- 
magnetism, when the system size increases, can be un- 
derstood with extreme statistics. This length scale has 
been confirmed with exact ground state calculations. The 
change of magnetization at the droplet excitation is nat- 
urally of "first-order" -kind. 

Above the break-up length scale we have studied the 
magnetization and susceptibility with respect to a con- 
stant external field. The behavior of the magnetization 
and the susceptibility is continuous and smooth and does 
not have any indications of a transition or a critical point, 
in agreement with the expectations of a continuously 
varying magnetization around H = 0, and a paramag- 
netic ground state. We thus conclude that the correct 
way of looking at the susceptibility is to study it with re- 
spect to the external field and above the break-up length 



scale instead of as a function of the random field strength, 
when the first order character of the break-up length scale 
may among others cause problems. 

However, we are able to find another critical phe- 
nomenon in the systems, in their geometry. For square 
lattices sites do not have a spanning property in ordinary 
percolation, when the occupation probability is one half. 
This corresponds to the random field case with high ran- 
dom field strength value without an external field. When 
an external field is applied and the random field strength 
decreased, a percolation transition can be seen. The tran- 
sition is shown to be in the standard 2D short-range cor- 
related percolation universality class, when studied as a 
function of the external field. Hence, the correlations in 
the two-dimensional random field Ising magnets arc only 
of finite size. We also want to point out that in these kind 
of systems, the random field strength is a poor control pa- 
rameter and the systems should be studied with respect 
to the external field, and after that map to the random 
field strength. By doing so we have been able to find a 
critical random field strength value, below which the sys- 
tems are always spanning even without an external field. 
When the percolation transition is studied without an 
external field and tuning the random field strength lots 
difficulties are encountered. This might cause puzzling 
consequences when studying the character of the ground 
states, not only because of the bad control parameter, 
but also because the percolation correlation length may 
be mistaken for something as the magnetization correla- 
tion length. Also note that the "true behavior" is seen 
only for system sizes large enough (L > Lb). 

The percolation character of the ground state struc- 
ture can be measured by the standard percolation frac- 
tal dimensional scaling for the mass of the spanning do- 
main. The existence of such a large cluster is not against 
the paramagnetic structure of the ground state, since the 
fractal dimension is below the Euclidean dimension. In 
order to be consistent with the Aizcnman-Wchr argument 
in the zero-external field limit the spins in the opposite 
direction from the external field may form the spanning 
cluster at low random field strength values. In fact we 
have found cases of finite systems for H — 0, where the 
magnetization of the system is opposite to the orienta- 
tion of spins in the spanning cluster. Notice that this 
does not imply that the critical lines i?c(A) actually cross 
each other at Ac continuing on the opposite side of the 
H = 0-axis (see Fig. 0) . By considering the red clusters it 
seems that in the TD-limit the spanning cluster should 
be broken up at = e, e — > since the field needed 
to fiip such a critical droplet should go to zero with L. 
Also, since the sum of the fields in the spanning cluster 
is shown to scale with the same fractal dimension as the 
mass, we conclude that the Imry-Ma argument does not 
work any more after the system has broken up on several 
domains. It works only for the first domain to appear. 

We have also generalized the red sites of the stan- 
dard percolation to red clusters in the percolation studied 
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here. A red cluster results from the energy minimization 
by flipping a whole cluster although only a single spin 
has been forced to be flipped, and breaks up the span- 
ning character of a percolating cluster. Actually the fi- 
nite size of the red clusters indicate also the presence of 
only short-range correlations in the systems. Such a lack 
of long-range correlations maybe explains, why we can 
see an "accidental" percolation phenomenon in a zero- 
temperature magnet whose physics is governed by the 
disorder configuration. The normal percolation univer- 
sality class is tightly connected to conformal invariance, 
which is most often destroyed by long-range correlations 
or randomness p^ ]. 

To finish the paper, we would like to raise some 
open questions related to the percolation behavior of the 
ground states of the two-dimensional random field Ising 
magnets. As noted, an interesting problem is the ex- 
act relation of the RFIM percolation to conformal in- 
variance. The percolation characteristics of the ground 
state might be experimentally measurable since the over- 
lap of the ground state and finite temperature magneti- 
zation should be close to unity for small enough temper- 
atures. The structure and relaxation of diluted antiferro- 
magnets in low external fields are suitable candi- 

dates: there one would presume it to be of relevance that 
there are large-scale structures present in the equilibrium 
state. In particular in coarsening, it is unclear how the 
eventual hull exponent of 4/3 would affect the dynamics. 
It would be interesting to see what kind of phenomena 
can be seen in the structure on triangular lattices since 
here Pc — 0.5 even in the ordinary site-percolation. One 
open question or application is the 3D RFIM. The perco- 
lation transition of the minority spins is expected to take 
place along a line in the {H, A) -phase diagram as well, 
since Pc — 0.312 for site-percolation in the case of the cu- 
bic systems most often studied numerically. Thus in low 
fields only one of spin orientations percolates whereas at 
high fields both do, see a review of 3D RFIM experiments 
in l^lj . The role of this transition is unclear also when it 
comes to the ferro- to paramagnet phase boundary, and 
the nature of the phase transition. 
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FIG. 1. An example of the ground state after the first exci- 
tation, — 200^, Gaussian disorder, A = 0.76. "Up" spins 
are drawn in white, and "down' spins in black. Note, that the 
system has periodic boundaries. 



Oh 
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FIG. 3. (color) Two examples of the ground state. "Up" 
spins in isolated domains are drawn in white, and "down" 
spins in black unless they belong to the spanning cluster 
(grey). The yellow spin is a seed of a red cluster, that breaks 
the spanning cluster. Periodic boundary conditions are used, 
and the spanning is checked in the vertical direction. The 
system size — 100^ and the random fields' standard devi- 
ations A = 3.0 and A = 1.2. 
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FIG. 2. The break-up length scale Lt versus inverse ran- 
dom field strength (1/A)'^ for bimodal and Gaussian disor- 
der (filled circles and empty squares, respectively) , calculated 
from Pfm{Li,) = 0.5. 
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FIG. 4. An example of the ground state, when an external 
field is applied. The system size L = 175 > Lb. Gaussian dis- 
order, A = 1.9. "Up" spins are drawn in white, and "down" 
spins in black. The external field H =0.0 (l.h.s. up), 0.1 
(r.h.s. up), 0.25 (l.h.s. down), and 0.5 (r.h.s. down). 
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FIG. 5. (a) The magnetization m versus the external field 
H for the system size = 175^ (all the tested system sizes 
L > Lb lie exactly on top of each other), and random field 
strength values A = 1.9, 2.0, 2.2, 2.6, 3.0, 3.5 and 4.5. Each 
point is a disorder-average over 5000 realizations and the er- 
ror bars are smaller than the symbols, (b) The external field 
values Hjni when the magnetization curves in (a) crosses the 
fixed magnetization values m — 0.25, 0.2, 0.1, -0.1, -0.2, -0.25, 
versus random field strength. The exp-lines are guides to the 
eye, and their prefactors are estimated using least-squares fit. 
(c) The data-collapse of (a) by scaling H with exp(— 6.5/A) 
estimated in (b). 
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FIG. 6. (a) The susceptibility, calculated as a fluctusr- 
tion of the magnetization, multiplied by the system size 
L^X ~ {fn^ — (ti)^) versus the external field H for A — 2.2. 
(b) The scaled susceptibility x/exp(7.3/A) versus the scaled 
external field H/ exp(— 6.5/A) for the system size = 175^ 
and random field strength values A = 1.9, 2.0, 2.2, 2.6, 3.0, 
3.5 and 4.5. Each point is a disorder- average over 5000 real- 
izations and the error bars are smaller than the symbols, (c) 
Same as (b) but only the positive external field values and in 
lin-log scale. The exp(— 0.2a;) line is to guide the eye. 
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FIG. 7. The phase diagram for the 2D RFIM with disorder 
strength A and an applied external field H. The 1/A = 
axis corresponds to the standard site-percolation, with the 
percolation occupation fraction pc = 0.593. The dashed lines 
define the percolation thresholds Hi;{l/A) for up and down 
spins, below and above which the systems are simple ferro- 
magnetic. The thick arrows denote two directions of which 
the percolation transition may be studied: the vertical one 
fixed A and varying H and the horizontal one H = and 
varying A. 
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FIG. 8. (a) The spanning probabilities of up spins n„p as 
a function of for A = 2.6 with I? € [20^ - 175^]. The 

data points are disorder averages over 5000 realizations the 
error bars being smaller than the symbols. The lines are 
sixth order polynomial fits, (b) The crossing points Hc{L) 
of the polynomials with the horizontal lines leads to the es- 
timate of the critical He by finite size scaling using L^^^" , 
V = 4/3. (c) The data-collapse with the corresponding criti- 
cal He = 0.00186 ± 0.0008. 
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FIG. 9. The scaled order parameter, probability of belong- 
ing to the up-spin spanning cluster, Poo/L-'^'", (3 = 5/36, 
V = 4/3 versus the scaled external field {H — He) / L~^^'^ , 
for A = 3.0 with e [20^ - 175^]. The data points are 
disorder averages over 5000 realizations the error bars be- 
ing smaller than the symbols. The corresponding critical 
Hc(A = 3.0) = 0.040 ± 0.001. 
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FIG. 10. (a) For each A the critical [//c(A)]^/* of up 
spin -spanning, where 5 — 2.05 ± 0.10. The data follows 
He ~ (A - Acf, where Ac = 1.65 ± 0.05. The details 
are as in Fig. bI Hc{A) spans almost two decades from 
He = 0.0028 ± 0.0008 for A = 2.0 to He = 0.1891 ± 0.002 
for A — 4.5. (b) The other scenario is shown, with 
He ~ A^ exp(-13/A2 - 4) (see text). 
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FIG. 11. (a) The spanning probabilities of up or down spins 
n as a function of A for J7 = with € [20^ - 1000^]. The 

data points arc disorder averages over 5000 realizations the 
error bars being smaller than the symbols. The lines arc tenth 
order polynomial fits, (b) The crossing points Ac{L) of the 
polynomials with the horizontal lines of 11 =0.1, 0.15, 0.2, 
0.25, 0.3, 0.4, and 0.5 versus the scaled system size. The es- 
timates of the critical Ac are searched by finite size scaling 
using Ao(L) = Ac(H-CiL-''''')(l + C2L"'/''=^) and the data 
is plotted as Ac(L) versus C\L~^^" . (c) The critical random 
field values Ac with respect the 11 values they are estimated 
from. The corresponding correlation length exponents l/i>, 
which are used in (b), are shown, too. (d) Another scenario: 
The critical random field values Ac with respect to system 
size. The lines are least-squares fits of form L ~ exp(C/A"), 
where C is a free parameter, for different 11 values from (a). 
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FIG. 12. The average mass of spanning clusters for bimodal 

field randomness A = 25/13 up to system size L = 470. The 
plot shows also the sum of the random fields of the sites be- 
longing to the same clusters. The 2d percolation fractal di- 
mension D =91/48 is indicated with a line. 
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FIG. 13. (a) The average number of red clusters, (Nnc), 
as a function of the system size, = 20^ — 200^, for A = 1.3, 
1.4, 1.5, and 1.6 with H = 0. Here the number of roahzations 
is N — 200. The smaUer the A the larger the amplitude of 
(Nrc), since the sizes of red clusters become larger. The line 
L^i"^ is a guide to the eye. (b) The masses of the red clusters 
(Mrc) with respect to the system size L. (Mrc) does not 
depend on the system size as seen in the figure. 
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